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We report the creation of ultralong-range Sr 2 molecules comprising one ground-state 5s 2 1 So 
atom and one atom in a 5sns 3 Si Rydberg state for n ranging from 29 to 36. Molecules are created 
in a trapped ultracold atomic gas using two-photon excitation near resonant with the 5s5p 3 Pi 
intermediate state, and their formation is detected through ground-state atom loss from the trap. 

The observed molecular binding energies are fit with the aid of first-order perturbation theory that 
utilizes a Fermi pseudopotential with effective s-wave and p -wave scattering lengths to describe the 
interaction between an excited Rydberg electron and a ground-state Sr atom. 


Low-energy scattering of the nearly free, excited elec¬ 
tron in a Rydberg atom from a ground-state atom can 
bind the two atoms into an ultralong-range Rydberg 
molecule 00 . The resulting internuclear spacing is on 
the order of the size of the Rydberg atom, which scales 
with the principal quantum number n, as n 2 and can 
exceed one micrometer. This class of molecules has at¬ 
tracted significant attention because it demonstrates a 
new mechanism for chemical bonding and the molecules 
possess surprising features, such as the presence of a per¬ 
manent electric dipole moment, even in the homo nuclear 
case Q. Here, we report the creation and theoretical 
description of ultralong-range Sr 2 molecules involving a 
bsns 3 Si Rydberg atom. 

Ultralong-range molecules were originally predicted 
theoretically [l| and were subsequently observed in Rb 
Q and Cs Q. The original observations were of dimers 
involving spherically symmetric S Rydberg states [2], but 
now measurements have been extended to anisotropic P 
[Bj] and D [b] Rydberg states and to molecules compris¬ 
ing one Rydberg atom and as many as five ground state 
atoms 0| . 

There has been increasing interest in ultracold Ryd¬ 
berg gases of alkaline-earth metal atoms because of sev¬ 
eral new possibilities introduced by their divalent elec¬ 
tronic structure. The principal transition of the Ryd¬ 
berg core is typically in the visible range and can be 
used to drive auto-ionizing transitions (§j, to image Ryd¬ 
berg atoms or ions [ 9 }, and to provide oscillator strength 
for magic-wavelength optical trapping of Rydberg atoms 
0- Doubly excited states serve as strong perturbers 
of Rydberg states, creating a much richer assortment 
of electronic configurations than found in alkali-metal 
atoms. The existence of triplet and singlet excited levels 
provides many Rydberg series, giving access to attractive 
and repulsive interactions Two-photon excitation to 
triplet Rydberg levels via the intermediate 3 Pi state, as 
used here, can also reduce the overall decoherence from 
light scattering for a given strength of optical coupling to 
the Rydberg level as compared to two-photon transitions 
available in alkali-metal atoms 0. The results reported 
here represent the first experiments involving ultracold 


Rydberg atoms excited via intermediate triplet excited 
states. 

Within the framework of a two-active-electron approx¬ 
imation, one of the two valence electrons can be excited 
to a Rydberg state. The interaction between the excited 
Rydberg electron and a neighboring ground-state atom 
can be described using the Fermi pseudopotential BUI 


Tr , ^ A 2nh 2 AJkCR)} r . 
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where and R specify the positions of the Rydberg- 
atom valence electron and ground-state atom, respec¬ 
tively, measured from the Rydberg core. The momen¬ 
tum dependent s-wave and p-wave scattering lengths 
are A s (k ) and A p (k). The Rydberg-electron mo¬ 
mentum in a semiclassical approximation is hk(r) = 
yj2m e (e 2 /(47T£or) — E\f), where E & is the binding energy 
of the unperturbed Rydberg atom. This approximation 
is justified for highly excited 3 Si Rydberg states because 
the electron-electron interaction is relatively small. For 
the case of both weak molecular binding compared to 
the binding energy of the Rydberg electron and a molec¬ 
ular size much larger than the characteristic size of the 
ground-state atom |I], a first-order perturbative approx¬ 
imation of the interaction energy leads to the mean-field 
molecular potential 


V(R) = 


2irh 2 A s [fc(R)] 


J dr 2 |^(ri = R,r 2 )| 2 (2) 
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where |\h(ri, r 2 )| 2 is the valence-electron probability den¬ 
sity for the Rydberg atom. The small correction due 
to p-wave scattering has been dropped in Eq. O For 
A s < 0, the interaction is attractive and can bind ground- 
state atoms to a Rydberg atom in potential wells created 
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by the antinodes of the electron wavefunction. For Sr, 
the theoretically predicted value of the e _ -Sr scattering 
length is A s (k = 0) = —18 ao flEj . 

The creation of ultralong-range molecules requires ul¬ 
tracold temperatures so that thermal energies are lower 
than their small binding energies 10MHz). Also, 
high-density samples are necessary to ensure a sizeable 
probability of finding two atoms with separations less 
than the radial extent of a Rydberg electronic wave- 
function. We obtain these conditions using 84 Sr atoms 
confined in an optical dipole trap (ODT). This isotope 
has collisional properties favorable for evaporative cool¬ 
ing and the creation of high phase-space-density samples 
[16| . Details of the cooling and trapping are given in 
Ref. [13 • 

At the start of the excitation time, the atoms are held 
in an ODT formed by crossed 1064-nm laser beams with 
waists of 300 fi m x 440 /im x 38 /im resulting in an approx¬ 
imately oblate spheroidal trap with axial frequencies of 
11 Hz and 12 Hz, and a radial frequency of 158 Hz. Typ¬ 
ically 7 x 10 5 atoms are trapped at a temperature of 
200 nK yielding a peak density of p = 2.7 x 10 12 cm -3 . 
This corresponds to an average interparticle spacing of 
1/p 1 / 3 = 0.7 /im, which is about 10x the radius of the 
outer lobe of the Rydberg wavefunction 2n* 2 ao — 75 nm 
for n — 30 and n* = n — 5. We determine the bsns 3 S\ 
state quantum defect S = 3.372 ± 0.001 by fitting ob¬ 
served Rydberg lines between n = 24 to 36, and this 
value agrees well with Ref. 0. 

Atoms are promoted to Rydberg states through 
two-photon 5s 2 1 So-5s5p 3 Pi-bsns 3 Si excitation. The 
689 nm laser for the first step is detuned 170 MHz to the 
blue of the intermediate state to avoid scattering from 
the atomic line and associated molecular resonances. Ry¬ 
dberg states with n = 29 — 36 are reached with photons 
at 319 nm generated by frequency doubling the red out¬ 
put of a fiber-based optical parametric oscillator laser. 
Approximately 200 mW of UV power is available. The 
intensities of the red and UV light on the atoms are 
2.2 x 10 3 W/m 2 and 2.3 x 10 5 W/m 2 respectively. The 
UV and 689 nm lasers co-propagate with orthogonal lin¬ 
ear polarizations. This configuration excites a superpo¬ 
sition of m = +1 and m = — 1 3 Si Rydberg states. 

The frequency of the UV light is controlled by locking 
the 638 nm fundamental to an optical cavity stabilized to 
the 689nm laser, which is locked to the 5s 2 1 So~5sbp 3 Pi 
atomic transition. The UV frequency is scanned using an 
acousto-optic modulator in a double-pass configuration 
in the path of the 638 nm light en route to the stabiliza¬ 
tion cavity. The excitation time is precisely controlled 
using an acousto-optic modulator on the 689 nm beam, 
while the UV light is controlled by a slower mechanical 
shutter. After excitation, the atoms are released from the 
trap and the ground-state atom population is measured 
with time-of-flight absorption imaging on the 5s 2 1 5'o- 
bsbp X P\ transition at 461 nm. Excitation to atomic Ryd¬ 
berg or molecular levels is detected as ground-state atom- 
loss. The exposure time is adjusted for about 50% peak 


loss, which is ^ 10ms for atomic resonances and ^ 2s 
for molecular resonances. The ODT is left on during ex¬ 
citation, and we assume the AC Stark shift is the same 
for the atomic and molecular transitions in our quoted 
results. 

Typical atom-loss spectra are shown in Fig. [lj The 
spectra are relatively simple because of the closed-shell 
1 S'o electronic-ground state and lack of nuclear spin 
for 84 Sr. Atom loss is an indirect method of detect¬ 
ing excitation as compared to the traditional technique 
of pulsed-field ionization and charged-particle detection 
[UlIHzl, but it can still yield a high signal-to-noise ra¬ 
tio. A comprehensive understanding of the decay chan¬ 
nels of ultralong-range molecules is still lacking, but it 
is clear that a large fraction of excitations should lead 
to measurable ground-state atom loss. Fluorescent de¬ 
cay after either atomic or molecular excitation has a siz¬ 
able probability, approximately two-thirds, of creating 
long-lived 3 Po or 3 P 2 atoms. Such atoms may remain 
trapped but they are invisible to absorption imaging. 
The recoil energies of h 2 /2m^\\^kB = 0.24 pK and 
h 2 /2ra84A 2 20 VB = 1.1/iK are smaller than the 1.9 pK 
depth of the ODT. However, collisional processes in¬ 
volving ground-state atoms and tunnelling to small in- 
ternuclear separation can reduce the molecular lifetime 
significantly compared to the atomic lifetime for con¬ 
ditions similar to ours [l8|. Such tunnelling should re¬ 
lease enough energy to eject one or both atoms from the 
trap. For atomic excitation, it is possible that the den¬ 
sity of Rydberg atoms becomes high enough for inelas¬ 
tic Rydberg-Rydberg collisions to lead to atom loss 0- 
The density of molecules, however, is low enough that 
such processes are negligible. 

The threshold for the molecular-binding energy for 
each principal quantum number is determined by mea¬ 
suring the resonance position for atomic excitation to 
the bsns 3 Si state. This is done with a 10 ms excitation 
time to avoid saturating the transition, and it results in 
a ~ 800 kHz FWHM linewidth, which is likely limited by 
the UV laser linewidth. For such short excitation times, 
no molecular transitions are visible. For excitation times 
on the order of 1 s, however, clearly resolved resonances 
corresponding to molecular bound states appear to the 
red of the (highly-saturated) atomic line. No transitions 
are observed to the blue or further to the red of the re¬ 
gions shown. Data are fit to the wing of a Lorentzian 
to describe the atomic background plus a gaussian for 
each molecular line. Molecular-binding energies are de¬ 
termined by the frequency difference between molecular 
and atomic lines. 

Typical molecular linewidths are 800 kHz FWHM, 
again limited by the UV laser linewidth. Several spectra 
were recorded for each principal quantum number, al¬ 
ternating between measurement of atomic and molecular 
lines. The uncertainties in the molecular line positions 
with respect to the atomic line, ±150 kHz, are the sta¬ 
tistical la uncertainties in the mean for each group of 
measurements. DC Stark shifts of the transitions were 
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FIG. 1: (Left) Atom-loss spectra and (Right) calculated potentials and wavefunctions,Rx^(R), for n=30 (top), 33 
(middle), and 36 (bottom). The bars at the top of the molecular spectra indicate the positions of the 
theoretically-predicted binding energies of states bound by > 1 MHz and the curves at the bottom indicate their 
calculated excitation strength. The origin of each frequency axis is set to the center of the atomic excitation 
spectrum (see text). 


measured for n = 36 by applying DC electric fields of up 
to about 0.5 V/cm. The extracted atomic and molecular 
DC polarizabilities were equal at our level of precision 
±0.5 MHz/(V/cm ) 2 and calculated theoretically to be - 
4.5 MHz/(V/cm) 2 . 

For the present range of n and our relatively low den¬ 
sities, the production rate for trimers is expected to be 
much lower than for dimers umi and therefore difficult 
to detect with our current methods. The most deeply 
bound level observed for each principal quantum number 
is assigned to the vibrational ground state of one atom 
in the potential well formed by the outermost lobe of 
the Rydberg electron wavefunction. This is confirmed by 
calculations of the potentials and molecular wavefunc- 
tions (Fig. [ljRight)), which show that this state is well 
localized in this well. The more weakly bound levels cor¬ 
respond to vibrationally excited states, which are delo¬ 
calized across several lobes of the electron probability 
density. 

The electron wave function T(ri,r 2 ) of the 3 S\ Ryd¬ 
berg atom can be calculated by numerically diagonalizing 
the Hamiltonian within the two-active-electron model. 
Spin-orbit interaction is included, but its effects are small 
|2Q| . The model potential is fitted to reproduce the mea¬ 
sured energy levels in the singlet sector [ 22 ], and yields 
the quantum defect 5 = 3.376 for 3 Si states, which agrees 
well with the measured value. The calculated wave func¬ 


tions for 3 Si Rydberg states with n ~ 30 are dominated 
by a single configuration 


T(ri,r 2 ) ~ (05s(ri)^ns(r 2 ) - V±s( r i)</>5 S ( 1 * 2 )) (3) 


where 05 S (r) and ip ns (r) are the wave functions of the 
5s state for Sr + ions and the Rydberg ns state for Sr 
atoms, respectively. The contributions from the other 
configurations are smaller than 0.01%. Considering that 
the wave function of the 5s state is rather localized (i.e. 
0 5 s (r) ~ 0 for r > 20 a.u.), the molecular potential (Eq. 
El around R ~ 1000 a.u. can be evaluated similarly to 
that for single-electron systems, 


v(R) = 


6tt h 2 Ap(k) 


m e 


m e 


|VV± S (R)| 2 . 


(4) 

By solving the Schrodinger equation associated with 
the molecular potential, the binding energies and the 
wave functions for Rydberg molecules are obtained. The 
effective scattering lengths are taken to be 


A s (k) = A s (k = 0 ) ± ^ ak 2 A p (k) = A p (k = 0 ). (5) 
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The measured value, a = 186ag, is used for the polar¬ 
izability [23| . It is known 0 that a non-perturbative 
Green’s function calculation correctly reproduces the 
measured molecular energy levels using the true zero- 
energy s-and p-wave scattering lengths. Within the 
first-order perturbative approximation, however, good 
agreement can be obtained by considering the scatter¬ 
ing lengths as effective fitting parameters. For stron¬ 
tium, while the calculated s-wave scattering length is 
A s (k = 0) = —18 a 0 0 , an effective scattering length of 
A s (k = 0) = —13.2ao for n = 30 to —13.3ao for n = 36 
yields good agreement with the measured energy levels, 
especially the most deeply bound molecular states. Al¬ 
though the depth of the deepest well in the molecular 
potential linearly scales with A s (k = 0) (Eq. [4]), the 
additional contribution from p-wave scattering becomes 
non-negligible around the nodes (i.e. |'0ns(R)| — 0) and 
affects the energies of the weakly bound levels and their 
density of states. They are optimally fitted using the 
value of A p (k = 0) 8.4ao- 

The molecular formation rate can be calculated as 
T v oc |(T, x^lTjTo, Xo)| 2 where T is the transition ma¬ 
trix for two-photon absorption, and To, Xo, Xv are the 
ground-state-electron wavefunction, the pair distribution 
of ground-state atoms, and the z^-th vibrational wave- 
function of the Rydberg molecule, respectively. Within 
first-order perturbation theory, the electronic part of the 
transition matrix is independent of the molecular states. 
Therefore, the excitation rate is reduced to the Franck- 
Condon factor 


oc 


/ 


dRR 2 Xu{R)xo{R) 


(6) 


Moreover the pair distribution xo(-R) is approximated as 
constant at large distances and can also be factored out. 
The excitation rate T v is convolved with a gaussian distri¬ 
bution and the obtained excitation spectra are included 
in Fig.ffl The agreement between the measured and cal¬ 
culated spectra is very good. The calculations reveal the 
existence of more molecular levels, but their excitation 
strengths are too weak to detect in the measurements. 

Figure [2] shows the n scaling of the observed binding 
energies together with the calculated values. The binding 
energies of deeply bound states follow the approximate 
— 1 /(n — S) 6 scaling seen in Rb [7]. Gaps in laser coverage 
prevented measurement of spectra for n = 31 and n = 
32. Deviations from the scaling are evident for weakly 
bound states and are clearly seen in high-quality spectra 
for n = 30. The energy scaling reflects the scaling of 
the probability density |'0ns(lT)| 2 of the Rydberg electron 
at the location of the unperturbed Sr atom. The wave 
functions of the most deeply bound states for 29 < n < 
36 are confined in the most outer well of the molecular 
potential (Fig. [TJ. Therefore, the corresponding energies 
scale with the depth of the well by — 1 /(n — S) 6 as the 
quantum number n varies. When the first-excited state 
of this potential well is nearly degenerate with the lowest 


Principal Quantum Number, n 

29 30 31 32 33 34 35 36 



FIG. 2: Scaling of observed molecular binding energies, 
showing — 1 /(n — S) 6 scaling for higher quantum 
numbers and more deeply bound levels. 


energy state of the adjacent well, the coupling between 
two states yields an energy splitting, A. If a molecular 
state confined in a single well for a given n experiences 
tunneling to an adjacent well for another value of n, the 
molecular energy does not scale with the potential depth 
showing a deviation ^ A from the —1 /(n — S) 6 scaling. 
Such a transition in the wave function is seen for some 
excited vibrational levels (Fig. [TJ). For example, although 
the first excited vibrational state’s wave function is rather 
confined within the first two outer wells for the given 
range of n, the second-excited state tunnels ever deeper 
into inner wells as n increases. Thus the deviation from 
the scaling is prominent for the latter (Fig. [TJ). Note that 
data obtained for n = 29 were of inferior quality, but we 
were able to identify the three lines indicated in Fig. [2] 
We have presented the first observation of ultralong- 
range Rydberg molecules in Sr 2 formed by photoexcita¬ 
tion to the bsns 3 Si Rydberg state for 29 < n < 36. The 
observed lines are well described using a Fermi pseudopo¬ 
tential approach to calculate perturbative molecular po¬ 
tentials and yield the effective s -wave and p-wave e - -Sr 
scattering lengths A s = — 13.2ao and A p = 8.4ao- 

This work represents the first study of ultralong-range 
Rydberg molecules in a divalent-atomic system, and it 
opens new directions in this emerging research area. Dou¬ 
bly excited electronic states give rise to dramatic pertur¬ 
bations of the Rydberg states in divalent atoms. This 
should lead to new types of ultralong-range molecules 
with mixed electronic character that arise from degen¬ 
eracies of pairs of Rydberg levels of different angular 
momenta. This may also lead to very strong transition 
strengths for production of molecules with large dipole 
moments. If one can form Rydberg molecules with high 
electronic angular momenta, it might be possible to opti¬ 
cally trap them using the oscillator strength of the ionic 
core 0. High densities of atoms in metastable triplet 
levels can be created in these systems [24|, allowing the 
formation of Rydberg molecules in which triplet atoms 
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serve as the “ground-state” atoms. Spectroscopy of these 
molecules will probe the low-energy scattering of elec¬ 
trons from the metastable states including measuring the 
electron-triplet scattering length, which should be sensi¬ 
tive to the greater polarizability compared to closed-shell, 
ground-state atoms. 

Detection of ultralong-range molecules with atom loss, 
as demonstrated here, greatly simplifies the required ex¬ 
perimental apparatus compared to charged particle de¬ 
tection, which may open the way for study of these ex¬ 
otic molecules in many other species of atomic gases. 
The measurements reported here also represent an im¬ 


portant step towards future experiments with interacting 
alkaline-earth Rydberg atoms because molecular excita¬ 
tions represent loss channels that need to be avoided. 
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